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We consider a huge quantum system that is subject to the charge superselection rule, which 
requires that any pure state must be an eigenstate ol the total charge. We regard some parts of 
the system as "subsystems" Si, S2, • • •, Sm, and the rest as an environment E. We assume that one 
does not measure anything of E, i.e., one is only interested in observables of the joint subsystem S 
= Si + S2 + ■ • ■ + Sm. We show that there exist states |$)tot with the following properties: (i) The 
reduced density operator TrE (|$)tot tot{4'|) is completely equivalent to a vector state |(/j)s s{'p\ of 
S, for any gauge-invariant observable of S. (ii) \'fi)s is a simple product of vector states of individual 
subsystems; \ip)s = \C^^^)i\C'^'^'^)2 ■ ■ ■ , where IC'*^)* is a vector state in Sfe which is not an eigenstate 
of the charge in Sfc. Furthermore, one can associate to each subsystem Sfc the vector state |C''°-')fc 
and observables which are not necessarily gauge invariant in each subsystem, and \C''''^)k is then a 
pure state. These results justify taking (a) superpositions of states with different charges, and (b) 
non-gauge-invariant operators, such as the order parameter of the breaking of the gauge symmetry, 
as observables, for subsystems. 

PACS numbers: 03.65.Bz, ll.15.Ex, 03.75.Fi, 74.20.-z 

In quantum theory, some superpositions of states are not permitted as pure states [Q. In particular, the charge 
superselection rule (CSSR) forbids coherent superpositions of states with different charges Q. Namely, any pure 
state must be an eigenstate of the total charge iVtot. However, it is customary to take such superpositions when one 
discusses the breaking of a gauge symmetry. Superconductors and superfluids are typical examples. If the system size 
V is infinite, this does not conflict with the CSSR because Atot becomes ill-defined as 1/ — > 00, and the CSSR becomes 
inapplicable. In real physical systems, however, phase transitions practically occur for finite {V < +00) systems as 
well. In particular, phase transitions have been observed in relatively small systems, including small superconductors 
1^, Helium atoms in a micro bubble |^], and laser-trapped atoms 0|. The meaning of the symmetry breaking in 
such systems has been a subject of active research j|-|8|. The purpose of this paper is to present a general discussion 
which justifies taking coherent superpositions of states with different charges for finite quantum systems subject to 
the CSSR, such as charged particles and massive bosons. Furthermore, we also justify taking non-gauge-invariant 
operators such as the order parameter of the breaking of the gauge symmetry, as observables of subsystems. 

We consider a huge quantum system that is subject to the CSSR. We regard some parts of the system as "subsys- 
tems" Si, S2, ■ ■ •, Sm, and the rest as the environment E. We assume the usual situation where (i) E is much larger 
than the joint subsystem S = Si -I- S2 -I- • • • -I- Sm, and (ii) one is not interested in (thus, one will not measure) degrees 
of freedom of E, i.e, one is only interested in S (or some parts of S). The Hilbert space Titot of the total system is the 
tensor product of the Hilbert spaces of Si, S2, • • •, Sm and E; Htot — Hs He, where Hs = Hi H2 ■ ■ • Ti-M ■ 
The total charge (in some unit) A^tot is the sum of the charges of Si, S2, • • •, Sm and E; 

M 

Atot = ^iVfe + 7VE. (1) 
fe=l 

Products of eigenfunctions |A^ini)i, • • •, {NmUm) m , and |A^e^)e of Ai, • • •, Nm, and A^e, respectively, form complete 
basis sets of Htot- Here, Nk (k — 1, 2, • • • , M) and Ae are eigenvalues of Afe and A^e, respectively, and and £ denote 
additional quantum numbers. 

The CSSR requires that any pure state of the total system must be an eigenstate of A'tot, i-e., superposition is 
allowed only among states with a fixed value of the eigenvalue A'tot of Atot . Consider the following state that satisfies 
this requirement: 

l*>tot - E • • • E E ^nL ■ ■ ■ ^^mIm^S • • • \NMnM)M\Ntot ~ Ns,£)e, (2) 

where A's — ^k, and the superposition coefficients are normalized as 

Ei^^'«i' = Ei<'.i' = i- (3) 
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For this state, the probabihty of finding iVtot — bosons in E takes almost the same vahies for all Ng such that 
\Ns — {Ns)\ < (JA^g)^/^ This property seems natural for a huge environment. Since we assume that one will not 
measure degrees of freedom of E, we are interested in the reduced density operator ps of S, which is evaluated as 



Ps = TrE(|$)tot tot($|) 

Ev^ f ^(1) MM) ^(M)* ^(1)* 

N[,n[ Afjj.n'jj- -/Vi,ni NM-nM 

X \N[n'^)i--- I A^m^i'm) Af M {Nmum | • • • i (iVmi | . (4) 
We can easily show that (ps)^ ^ ps unless 



for some set of numbers ■ ■ ■ , Nq^'^K We exclude this trivial case (where the CSSR is satisfied in each subsystem) 
from our consideration. Then, (ps)^ 7^ ps, and one may say that pg represents a mixed state. However, the relation 
(ps)^ 7^ Ps only ensures that for any vector state \(p)s (G Hs) there exists some operator Ss (on Hs) for which 

Trs (^psSs) ^ s{v\^s\^)s- (6) 

Note that such a general operator Ss is not necessarily gauge-invariant. Hence, Sg might not be an observable, which 
must be gauge- invariant. In fact, we first show that ps is equivalent to a vector state |(p)s for all gauge-invariant (thus 
physical) observables on Tig. This statement might not sound surprising because a vector state is not necessarily a 
pure state [Here, we use the precise definition of pure and mixed states rather than misleading definitions 

such as p^ = p and p^ ^ p.] In fact, the equivalence of |(p)s to ps means that the vector state \ip)s is a mixed state. 
In other words, \f)s is a vector state in a reducible representation of the algebra of gauge-invariant observables [ p^ . 
Actually, we first show a stronger statement: \ip)s is a simple product of vector states of individual subsystems; 

Ps is equivalent to |(p)s = \C'--^^)i\C^^'')2 ■ ■ ■ \C^^^^)m for any gauge- invariant observables in Hs, (7) 
where \C'^''^)k is a coherent superposition of states with different charges; 

\C^'''^)k=J2^Ni\Nn),. (8) 

N.n 

To see this, we recall that one will not measure degrees of freedom of E. This means that one measures only observables 
which take the following form; 

As » iE, (9) 

where ^s is an operator on Tis, and Ie denotes the unity operator on TYe- Note that ^s ® Ie must be gauge-invariant 
because of the gauge invariance of the total system, hence ^s must also be gauge- invariant. This requires that iVs 
(= J2k ^k) should be conserved by the operation of Ag. Hence, the matrix elements of As should take the following 
form; 

i(iVi7ii| • • • MiNMUMl As \K,n\,)M ■ ■ ■ \N[n[), - <5^,+...+jv.„wi4-...+Jv;,, A;^!::-;:::;^f (10) 
Hence, the expectation value of ^s for ps is evaluated as 
(As) = Trs (psis) 

\ " ^(1) ^(M) ^(M)* ^(1)* r. ANini---NMnM 

- 2^ "2^ 2^'" 2^ "^Wlni ■■■'-'w;,n^^A'MnM ■■■'"Afini°^^i + ---+JVM,W( + -..+Ar;,^^|„'^...Ar;^ 

= s(^|is|^)s, (11) 

and Eq. (0) is proved. The point is that Eqs. (^ and (^ contain the same Kronecker's delta, Jatj-i \-Nm,n[^ t-w^,^- 

Although this factor in Eq. (^) makes ps different from |(p)s s(<p|j the difference becomes totally irrelevant to {As) 
because of the same factor in Eq. ([l^). (Recall that {5ni^ \-Nm,n['\ hw^)^ = "^Wih ^Nm,n[^ ^n'^-) 
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Moreover, we can easily show that ps-Sm = Tr^f [ps] is also equivalent to a vector state, |C(^))i |C^^^)2 • ■ ■ IC^^^^W-i, 
where Tr^ denotes the trace operation over Hm- This result is easily generalized: For states of the form of (||), its 
reduced density operator for any set of subsystems is completely equivalent to a vector state, which is a simple product 
of vector states of the individual subsystems, if one measures only gauge- invariant observables of the subsystems. Note 
that there is no 'entanglement' between any subsystems In general, the absence of entanglement among sub- 

systems means the 'separability', i.e., it is possible to control quantum states of individual subsystems independently 
by local operations Hence, one can control the state \C^'''^)k of each subsystem by local operations, without 

perturbing the other subsystems |l4[ . 

For example, for interacting many bosons, which are confined in a large but finite box, one may take |C^*''')/£ to be 
the 'coherent state of interacting bosons' (CSIB), which is defined by |]8|, p^ , p^ 

oo N 

|a,,G), = e-l"'=l^/2^^|iV,G),, (12) 

where ak = |afc|e*'^'' is a complex amplitude, and | A^, G)k denotes the ground state that has exactly N bosons, which we 
call the 'number state of interacting bosons' (NSIB) One may also take jC*^'^); as \ai, G)i, where ai = \ai\e^'^' . 

Then, ps is equivalent to 

|^)s = ---\ak,G)k---\ai,G)i--- (13) 

This state has an almost definite value of the relative phase (j)ki between two condensates in Sfe and Sj. To see this, 
we may define the operator (acting on Hk <8) Hi) corresponding to e*"^*' by 

g#fci ^ gi4>k (e*"^')^. (14) 
Here, 6*"^*=, acting on Ti^, is not an exponential of some phase operator (j)k, but is defined by 

^>' = {bl,bok + ir'^^bok, (15) 

where bgk denotes the operator bo, which is defined in Ref. [l5|] as a nonlinear function of free operators, for S^. This 
operator is a linear combination of the cosine and sine operators of Rcf . . In the same way, e'"^' , acting on Hi , is 
defined using 69/ , which is 60 for S;. By similar calculations as in Ref. [|l5|, we can show that e*"^'', (e"^')^, and e^'t'^i 
have almost definite values, and we can regard 0fc; ~ 0fc — (f>i. Hence, the state (^3|) has the almost definite value of 
the relative phase. It should not be confused with states of the following form, which have been frequently discussed 
in the literature |l^ ]; 

Ms = E E ■ ■ • ^'"""^"Rn, \Nk,G)k ■ ■ ■ e'(^-^'')*'i?jv-w. \N - Nk,G)i • • • , (16) 

N Nk 

where Cn is a complex coefficient and i?'s are real ones. This state also has an almost definite relative phase if i?'s are 
appropriately taken. However, since Sfe and S; are strongly entangled in this state p^ , one cannot control the state 
of Sfe by local operations on Sfe, without perturbing the state of S/. In contrast, this is possible for the simple product 
state \<f)s- Moreover, we now show that \(f)s allows us to treat observables of individual subsystems separately. 

As exemplified by 6*'^'=' above, is generally a sum of products of operators of subsystems [T^ ; 

As = J2AkAi---, (17) 

where Ak (Ai) is an operator on Hk (Hi), and the sum is over combinations of AkAi ■ ■ ■. For the product state \if)s, 
the expectation value of Ag is simply evaluated from the expectation values of individual subsystems as 

s(^|is|^)s = E fc(C'('='|ifc|C«)fe i(C(')|iz|C«)z • • • . (18) 

From Eqs. (^) and (p^, we can consider each subsystem separately from the other subsystems. Namely, for subsystem 
Sfe, we can consider that its state is the vector state \C^'^^)k, and that Ak is one of its observables. When the 
expectation value of As for the set of subsystems is necessary, it can^ be evaluated from the expectation values of 
individual subsystems, as Eq. (111). Note that the gauge invariance of As does not require the gauge invariance of Ak 
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of each subsystem: It rather requires that only gauge- invariant combinations of A^Ai ■ ■ ■ should appear in the sum. 
[Namely, each N^^ can vary by the operation of Aj^, whereas (= Ni^) is conserved by the operation of Ai^Ai ■ ■ •.] 

For example, 6*"^*=' is gauge invariant, whereas neither e*"^*" nor e**^' is. Hence, by considering Ak as an observable 
of Sfe, we can include non-gauge invariant operators into the set of observables in Ti^,, with the restriction that 
among various combinations ol A^Ai ■ ■ ■ we must take only gauge-invariant ones as physical combinations: Results for 
other combinations should be discarded. This formulation gives correct results for all gauge-invariant (thus physical) 
combinations. The vector state |C'-'^^)fe then becomes a pure state of an irreducible representation of the algebra of 
observables of Sfe , because non-gauge invariant observables are now included in the algebra ||l^ . We have thus arrived 
at non-gauge invariant observables and a pure state which is a superposition of states with different charges, for each 
subsystem. 

The restriction that we must take only gauge-invariant ones among various combinations of AkAi ■ ■ ■ has the 
following physical meaning: If A}, is a non-gauge invariant observable of S^, then its value can only be defined relative 
to some reference observable Ai of some reference system S;. When \ip)'s, takes the form of Eq. (p^, for example, we 

can consider that Sfc is in a pure state ja^, G)k, and that the phase factor e^'^*' is one of observables on Tik- However, 

like the classical phase factor, the quantum phase factor e^'^'' can only be defined relative to some reference (such 

as (e*"^')^ of S;). In contrast to the position observable, which also requires a reference system but it can be any 
system composed of any material, the reference system of the phase factor should contain particles of the same kind 
as the particles in S^. Hence, the system of interest (S^) and the reference system (Sj) should be subsystems of a 
larger system of the same kind of particles. When the larger system contains a huge system(s), of which we are not 
interested in, we may call it an environment E, and the present model is applicable. When one is only interested in Sfc, 
the reference system S; may be considered as a part of the external systems, which include an apparatus with which 
one measures (jiu- In this case, the reference system S; can be regarded as a part of the measuring apparatus. (In 
the case of measurement of the position, for example, a material which defines the origin of the position coordinate 

can be considered as a part of the measuring apparatus of the position.) Although e*'^'' is not gauge invariant, results 

of measurement is gauge invariant because e*"^* (e*"^')^ is gauge invariant. In general terms, although results of any 
physical measurements must be gauge invariant, it does not necessarily mean that any observables of a subsystem 
must be gauge-invariant, because the gauge invariance of A^Ai (which ensures the gauge invariance of results of 
measurements) does not necessarily mean the gauge invariance of Ak . 

Finally, we present some significant applications of the present results. First, Our results give a natural interpre- 
tation of the order parameter O of a finite system, which can exhibit the breaking of the gauge symmetry in the 
infinite- volume limit. Although one might suspect that O must be excluded from observables since it is not gauge 
invariant, our results allow one to include it among observables as one of A^s. Second, one might also suspect that 
the CSSR would forbid a pure state which has a finite expectation value of O because its expectation value vanishes 
for any eigenstate of N. However, our results show that such a pure state is allowed as a state of a subsystem, |C''^))fc, 
which is not entangled with states of the other subsystems The absence of entanglement allows the preparation of such 
a pure state of a subsystem without perturbing the other subsystems (with perturbing the environment only). For 
interacting many bosons, for example, the order parameter of the condensation is usually taken as the field operator 
■0 of bosons. The present work justifies taking such a non-gauge-invariant operator as an observable of a subsystem. 
The expectation value of ij} is finite only for superpositions of states with different numbers of bosons. The present 
work justifies taking such a superposition as a pure state of a subsystem. More concretely, the expectation value of 

is finite for the CSIB, {a,G\ip\a,G) = 0(1), whereas it vanishes for the NSIB, (iV, G|^|7V, G) = §,0. Although 
both the NSIB and CSIB have the off-diagonal long-range order [plp^], the gauge symmetry is broken, in the sense 
that (■!/') = C(l), only for the CSIB Since the state vectors of the NSIB and CSIB are quite different, they have 
different properties. The most striking difference is that the NSIB decoheres much faster than the CSIB when bosons 
have a finite probability of leakage into a huge environment, whose boson density is zero Q. Moreover, it was shown 
that the CSIB has the 'cluster property,' which ensures that fluctuations of any intensive variables are negligible, in 
consistency with thermodynamics ||2^. Although the CSIB may look against the CSSR, the present work justifies 
taking it as a pure state of a subsystem, which, unlike the NSIB, is robust, symmetry breaking, and consistent with 
thermodynamics. In particular, generalizing Eq. (p^, one can take CSIB's for all subsystems. In this case, the joint 
subsystem S is in a pure state in which each subsystem is in a CSIB. Namely, the state vector of S behaves locally as a 
CSIB, which is robust, symmetry breaking, and consistent with thermodynamics. Moreover, unlike Eq. (p^), one can 
controlc states of individual susbsystems independently by local operations. These justify macroscopic theories such 
as the Ginzburg-Landau theory and the Gross-Pitaevskii theory, which assume, sometimes implicitly, that the order 
parameters can be defined locally, and that their fluctuations are negligible, and that the state vectors are robust 
against weak perturbation from environments, and that local operations do not cause global changes. We consider 
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that such a locally-CSIB state should be (a good approximation to) a quantum state of real physical systems at low 
temperature, except for extreme cases such as bosons in a perfectly-closed box at a ultra-low temperature. Finally, 
we remark that we have not assumed that the particles in the environment are in (or, not in) a condensed phase. 
We thank H. Tasaki and M. Ueda for discussions. 
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